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Dynamic  fracture  analysis  of  concrete 
structures  necessitates  a triaxlal  stress-strain 
relation  that  describes  gradual  strain-softening 
with  reduction  of  tensile  stress  to  zero.  A new 
model  which  does  that  and  is  applicable  under 
general  loading,  including  rotating  principal 
stress  directions,  is  proposed.  It  is  based  on 
accumulating  stress  relaxations  due  to  micro- 
cracking from  the  planes  of  all  orientation  within 
the  microstructure.  Comparisons  with  tensile  test 
data  are  given. 


Introduction 

Fracture  analysis  of  certain  brittle  heter- 
ogeneous materials,  such  ns  concretes  and  many 
rocks,  requires  consideration  of  progressive  micro- 
cracking in  the  fracture  process  zone  as  the  frac- 
ture is  being  formed.  This  type  of  fracture  may  be 
efficiently  modeled  with  the  crack  band  approach, 
in  which  the  material  behavior  in  the  fracture 
process  zone  is  described  by  a strain-softening 
triaxlal  stress-strain  relation,  provided  that  the 
strain-softening  behavior  is  associated  with  a 
zone  of  a certain  characteristic  width  that  is 
treated  as  a material  property  or  is  determined  in 
advance  by  stability  analysis.  A suitable  tri- 
axial  stress-strain  relation  of  the  total  strain 
type  (deformation  theory  type)  has  been  recently 
formulated  and  has  been  shown  to  lead  to  satis- 
factory agreement  with  essentially  all  existing 
fracture  test  data  available  in  the  literature 
[1,  2].  This  stress-strain  relation  is,  however, 
limited  to  situations  in  which  the  direction  of  the 
maximum  principal  stress  does  not  significantly 
rotate  during  the  fracture  formation.  This  is  not 
so  in  certain  important  situations,  especially 
various  dynamic  problems,  here,  a longitudinal 
wave  may  produce  only  a partial  tensile  fracture 
(i.e.,  distributed  microcracking)  and  the  fracture 
may  be  completed  subsequently  when  a shear  wave 
arrives,  causing  a principal  tensile  stress  in  a 
different  direction.  For  such  situations  of  pro- 
gressive fracturing,  it  is  necessary  to  develop  a 
triaxlal  strain-softening  stress-strain  relation 
which  is  path-dependent  and  is  formulated  Incre- 
mentally. A model  called  raicroplanc  model  is  de- 
veloped to  fill  this  need.  We  propose  here  a 
model  in  which  the  constitutive  properties  are 
characterized  by  a relation  between  the  stresses 


and  strains  acting  within  the  microstructurc  on 
planes  of  various  orientation,  called  the  micro- 
planes.  This  formulation  Involves  no  tensirial 
Invariance  restrictions.  The  restrictions  can 
then  be  satisfied  by  a suitable  combination  of 
planes  of  various  orientation.  E.g.,  in  the  case 
of  Isotropy,  each  orientation  must  be  equally 
frequent.  Thus,  one  circumvents  the  difficulty  of 
setting  up  a general  nonlinear  constitutive  equa- 
tion in  terms  of  proper  invariants. 

The  idea  of  defining  the  inelastic  behavior 
independently  on  planes  of  different  orientation 
within  the  material,  and  then  in  some  way  super- 
imposing the  inelastic  effects  from  all  planes, 
appeared  in  Taylor's  work  (3)  on  plasticity  of 
polycrystalline  metals.  Batdorf  and  Budianski  [6) 
formulated  the  slip  theory  of  plasticity,  in  which 
the  stresses  acting  on  various  planes  of  slip  are 
obtained  by  resolving  the  macroscopic  applied 
stress,  and  the  plastic  strains  (slips)  from  all 
planes  are  then  superimposed.  The  same  super- 
position of  inelastic  strains  was  used  in  the  so- 
called  multilaminate  models  Df  Zienkicwic z ct  al. 
15]  and  Pande  et  al  |6.  7]  and  in  many  works  on 
plasticity  of  polycrystals.  While  the  previous 
works  dealing  with  plasticity  of  polycrystals  L 3 , 

6,  8,  9,  10-14]  or  soils  [15,1b]  the  stresses  on 
various  microplanes  were  assumed  to  be  equal  to  the 
resolved  macroscopic  stress,  this  new  model  uses  a 
similar  assumption  for  part  of  the  total  strains. 

Fundamental  Hypotheses 

The  resultants  of  the  stresses  acting  on  the 

raicroplanes  over  unit  areas  of  the  macroscopic 

continuum  will  be  called  the  microstresses  s.., 

ij 

and  the  strains  of  the  macroscopic  continuum  accu- 
mulated from  the  deformations  on  the  microplanes 


will  be  called  the  microstrnins , e 


With  regard 


to  the  interaction  between  the  micro-and  macro- 
levels, one  may  Introduce  the  following  basic 
hypotheses . 

Hypotehsls  I.  - The  tensor  of  macroscopic 
strain,  r^j>  is  a sum  of  a purely  elastic  macro- 

strain  e“.  that  Is  unaffected  by  cracking,  and  an 
inelastic  macrostrain  e^j  which  reflects  the 
stress  relaxation  due  to  cracking,  i.e. , 
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attjggigfr 


c . . * c . . + e.  . 
Ij  ij  i.l 


pieces,  since  mi croc racking  is  chiefly  concentrated 

there . 


Hero,  lalin  lower  case  subscripts  refer  to  car- 
tesian coordinates  x (i  = 1,2,3). 


In  llypotehsis  111,  the  relaxation  of  shear 
microstresses  s;  caused  by  the  shear  and  normal 


Hypothesis  11.  - The  normal  microstrain  e^ 

which  governs  the  progressive  development  of  crack- 
ing on  a inicroplane  of  any  orientation  is  equal  to 

the  resolved  macroscopic  strain  tensor  e..  for  the 

ij 

same  plane , 1 . e . , 

e = n . n . e . . (2) 

n i J ij 

in  which  n.  - direction  cosines  of  the  unit  normal 

t 

n of  the  inicroplane  and  the  repeated  latln  lower 
case  subscripts  indicate  a summation  over  1,  2,  3. 

Hypothesis  111.  - The  stress  relaxation  due 
to  all  microcracks  normal  to  is  characterized  by 
assuming  that  the  microstress  s^  on  the  microplane 

of  any  orientation  is  a function  of  the  normal 
microstrain  on  the  some  plane,  i.e., 

sfi  • (2n/3)  F(en)  (3) 

The  factor  (2n/3)is  introduced  Just  for  convenience, 
as  it  will  later  cancel  out. 

The  last  hypothesis  is  similar  to  that  made 
for  shear  microstresses  and  microstrains  in  the 
slip  theory  of  plasticity.  Hypothesis  11  is  how- 
ever opposite.  There  are  three  reasons  for  hypo- 
thesis II. 

1.  Using  resolved  stresses  rather  than  resolved 
strains  on  the  microplanes  would  hardly  allow  des- 
cribing strain-softening,  since  In  this  case  there 
are  two  strains  corresponding  to  a.  given  stress  but 
only  one  stress  corresponding  to  a given  strain. 

2.  The  microstrain6  must  be  stable  when  the 
macrostrains  are  fixed.  It  has  been  experienced 
numerically  that,  In  the  case  of  strain-softening, 
the  model  becomes  unstable  if  resolved  stresses 
rather  than  strains  arc  used. 

3.  The  use  of  resolved  strains  rather  than  re- 
solved stresses  seems  to  reflect  the  microstructure 
of  a brittle  aggregate  material  more  realistically. 
The  use  of  resolved  stresses  is  reasonable  for  poly- 
crystalllne  metals  in  which  local  slips  scatter 
widely  while  the  stress  is  roughly  uniformly  dis- 
tributed throughout  the  microstructure.  By  con- 
trast, In  a brittle  aggregate  material  consisting 

of  hard  Inclusions  embedded  in  a weak  matrix,  the 
stresses  are  far  from  uniform,  having  sharp  ex- 
tremes at  the  locations  where  the  surfaces  of 
aggregate  pieces  are  nearest.  The  deformation  of 
the  thin  layer  of  matrix  between  two  aggregate 
pieces,  which  yields  the  major  contribution  to  in- 
elastic strain,  is  determined  chiefly  by  the 
relative  displacements  of  the  centroids  of  the  two 
aggregate  pieces,  which  roughly  correspond  to  the 
macroscopic  strain.  The  microplanes  may  be 
imagined  to  represent  the  thin  layers  of  matrix 
and  the  bond  planes  between  two  adjacent  aggregate 


lion  is  probably  quite  good  for  very  small  crack 
openings,  since  it  has  been  deduced  from  test  data 
on  shearing  of  cracks  in  concretes  that  no  rela- 
tive shear  displacements  on  the  rough  Interlocked 
cracks  is  possible  before  a certain  finite  crack 
opening  is  produced,  and  that  the  shear  stiffness 
of  the  cracks  decreases  rather  slowly  as  the  crack 
gradually  opens.  One  must  admit,  however,  that 
Kq.  3 (Hypothesis  111)  is  also  justified  by  Its 
simplicity.  It  would  be  much  more  complicated  to 
assume  a general  relation  between  the  normal  and 
shear  microstresses  and  microstrains  on  each  plane. 

Tangent iol  Stiffness  Matr ix 

The  virtual  work  of  stresses  per  unit  volume 
may  be  written,  according  to  Kq.  1,  as  AW 
= a..  Sc,,  - 0.,  5e‘?,  + o.,  Se,  Summing  the 

ij  1.1  ij  ij  ij  ij 


virtual  work  due  to  and  5e^,  uc  further  have 

AW  = o?.  Ae'].  + s,,  6e..,  in  which  s..  is  the 
ij  ij  ij  ij  13 

macrostress  tensor  resulting  from  on  all  planes, 

and  o'\  is  the  stress  tensor  corresponding  to  c ^ . 

Since  both  expressions  for  6W  must  hold  for  any 

and  any  we  must  have  » o'^ . = o^. 

Rquilibrium  conditions  may  be  expressed  by 
means  of  the  principle  of  virtual  work: 


^ ■ 3 * aij 


6e, . = 2 
i.l 


/ sn  Sen  f (n) 

Js  n n 


in  which  S represents  the  surface  of  a unit  hem- 
isphere, the  factor  (4v/3)is  due  to  integrating 
over  the  surface  of  a sphere  of  radius  1,  and  JS 
- sinAdOdi  (Fig.  lb).  Note  that  we  do  not  need  to 
Integrate  over  the  entire  surface  of  the  sphere, 

since  the  values  of  o or  e are  equal  at  any  two 
n n 

diametrically  opposite  points  cn  the  sphere. 
Function  f (n)  defines  the  relative  frequency  of  the 
planes  of  various  orientations  n,  contributing  to 
inelastic  stress  relaxation. 

Substituting  F.qs.  2-3  into  JScj . A,  we  get 

o..  5e, , = / F(e  ) n.  n,  Ae.,  f(n)  dS,  and  because 
ij  ij  Js  n i ) ij 

this  must  hold  for  any  Se^,  we  must  have 


2ti  fu/2 


•n-n. 


F(e  ) n.  n.  f(n)  sinb  d$  dG  (5) 
n i j 


Furthermore,  according  to  Eq.  2,  dF(en> 

« F'(e  )de  = F' (e  )n.  n de.  , and  thus  differ- 
n n n k m km 

entiation  of  Eq.  5 finally  yields 

do  = Dc  . de.  (6) 

ij  ijkm  km 
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in  which 

:'2  ti/Mi/2  y 

l)?.,  ~ I a...  F*  (e  ) f (n)  Kiivii  d<J»  dQ,  with 

ijkm  * J ijk.ni  n 

0 J0  (7) 

a . , . - n . n , n.  n 

Ljkm  ijkm 

Pijkm  the  tangent  stiffnesses  of  the 

microplane  system.  Note  that  the  sequence  of  sub- 
scripts of  I)*:.,  is  immaterial;  therefore,  there 
i jkin 

are  only  six  independent  values  of  incremental 
stiffnesses.  Eq.  7 applies  to  initially  aniso- 
tropic solids.  For  isotropic  solids,  we  may  sub- 
stitute f (n)  = 1 . 

The  mathematical  structure  of  the  urgent 
model  may  be  geometrically  visualized  with  the 
rheologic  model  in  Fig.  la. 


The  stress-strain  relation  for  the  micro- 

planes.  relating  o to  t , must  describe  cracking 
it  n 

.ill  Lhe  w;iy  to  complete  fracture,  at  which  o 


reduces  to  zero.  In  view  of  the  kinematics 
visualized  in  Fig.  lb,  it  is  clear  that  <«  as  a 

function  of  c must  first  rise,  then  reach  a max- 
imum, and  then  gradually  decline  to  zero.  We 
choose  the  final  zero  value  to  be  attained  asymp- 
totically, since  no  precise  information  exists  on 

the  final  strain  at  which  a =0,  and  since  a 

11 

smooth  curve  is  convenient  computationally.  The 
following  expressions  were  used  in  computations 
(19]  (Fig.  lc)i 


-<ke  P) 

for  e > C:  c = E c e n 

n n a n 

lor  e < 0:  o ■ E e 

n n n n 


(ID 


The  compliance  corresponding  to  the  addi- 
tional elastic  strain  a'!,  must  satisfy  isotropy 
conditions,  and  so 


'ijkm 


9K 


--  i,  . S.  +•—-(«..  5.  - t 6..  i,  ) (8) 

a i)  km  2ga  lm  3 (J  hm 


in  which  Ka  and  G'1  are  certain  bulk  and  shear 
moduli  which  cannot  be  less  than  the  actual  initial 
bulk  and  shear  moduli  K and  G.  For  fitting  of 
test  data,  it  was  assumed,  with  success,  that 
1/Ga  = 0. 


Recalling  Kq.  1 (and  Fig.  la),  we  may  now 
write  the  incremental  stress-strain  relation  as 


d°ij  = Dljkm  <leij>  Wlth  ‘"i.jkJ  " 

l<DC\k*+ 


(9) 


_y  Applying  Eq.  7 to  elastic  deformations  (with 
f(n)  “ 1),  one  finds  that  the  matrix  in  Kq.  7 
always  yields  Poisson's  ratio  V = 1/A.  This  is  bo 
cause  the  slips  on  all  microplaacs  are  neglected. 
Since  v ■ 1/A  is  not  quite  true  for  concrete,  the 
additional  elastic  strain  must  be  used  to  make  a 
correction.  l.et  us  now  determine  the  value  of  K 
needed  to  achieve  the  desired  Poisson's  ratio  v. 
let  superscripts  c and  a distinguish  between  the 

values  corresponding  to  D^jgm  311,1  Vot 

uniaxial  stress  we  have  = <>u/9Ka  + o^/Ec 
and  e,.,  “ o,,/9Ka  - vco,,/Ec  in  which  uc  = 1/A 

Li.  11  11 

and  Ec  * 2ttE  /5»  E - F*(0)  = Initial  normal 
n n 

stiffness  for  the  rolcroplane.  Since  = - vcjj* 
wo  must  have 


Ka  = 

t + V 

— EC 

(for  v < vC) 

J. 

9(vC  - 

v) 

§ 

This 

is,  of 

course. 

under  the  assumption  t' 

1/Ga 

“ 0. 

in  which  F.  , k,  and  p are  positive  constants:  k a 
7 n 

1.8  x 10  , p *=  2. 


The  integral  in  Kq.  7 has  to  be  evaluated 
numerically,  approximating  it  by  a finite  sum: 


i jkm 


w [a.  F* 

ft  ljkm 


(e)l 


(12) 


in  which  ci  refers  to  the  values  at  certain  numer- 
ical integration  points  on  a unit  hemisphere  (i.e., 
certain  directions),  and  are  the  weights 

associated  with  Lite  integration  points. 


Since  in  finite  element  programs  for  incre- 
mental loading  the  numerical  integration  needs 
to  be  carried  out  a great  number  of  times,  a very 
efficient  numerical  integration  formula  is  needed. 
For  the  slip  theory  of  plasticity,  the  integration 
was  performed  using  a rectangular  grid  in  the  8-$ 
plane.  This  formula  is,  however,  computationally 
inefficient  because  the  integration  points  are 
crowded  near  the  poles,  and  also  because  in  the 
0— 4»  plane  the  singularity  arising  from  the  pole 
takes  away  the  benefit  from  a use  of  higher-order 
integration  formula. 


Optimally,  the  integration  points  .should  be 
distributed  over  the  spherical  surface  as  uniformly 
as  possible.  A perfectly  uniform  subdivision  is 
obtained  when  the  microplanes  normal  to  the  o- 
di ructions  are  the  sides  of  a regular  polyhedron. 

A regular  polyhedron  with  the  most  sides  is  the 
icosahedron,  for  which  N = 10  (half  the  number  of 
sides).  Such  a numerical  integration  formula  was 
proposed  by  Albrecht  and  Collatz  (18). 


Numerical  experience  revealed,  however,  that 
10  points  are  not  enough  whei.  strain-softening 
takes  place;  it  was  found  that  the  strain-softening 
curves  calculated  for  uniaxial  tensile  stresses 
oriented  at  various  angles  with  regard  to  the  in- 
directions significantly  differ  from  each  other, 
even  though  within  the  strain-hardening  range  the 
differences  are  negligible.  Therefore,  more  than 
10  points  are  needed,  and  then  a perfectly  uniform 
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Table 


1 - Direction  Cosines  and  Weights 
0(h®) 


for 


2 x 21  Points 


(Orthogonn 1 ) 


a 

X 

ft 

X2 

a 

x3 

Cl 

V 

i 

1 

0 

0 

0.02652141274 

2 

0 

1 

0 

3 

0 

0 

1 

II 

4 

0. 7071067812 

0.70/1067812 

0 

0.01993014153 

5 

0.7071067812 

-0. 7071067812 

0 

II 

6 

0.7071067812 

0 

0. 7071067812 

II 

7 

0.7071067812 

0 

-0. 7071067812 

II 

8 

0 

0. 7071067812 

0.7071067812 

II 

9 

0 

0.7071067812 

-0. 7071067812 

II 

10 

0.3679072746 

0.38790/2746 

0.8360956240 

0.02S07L242/2 

11 

0.3879072746 

0.3879072746 

-0.3360956240 

II 

12 

0.3879072/46 

-0.3879072746 

0.8360956240 

13 

0.3879072746 

-0.3879072746 

-0.8360956240 

11 

14 

0.3879072746 

0.8360956240 

0. 3879072746 

II 

15 

0.3879072746 

0.8360956240 

-0. 3879072746 

" 

16 

0.3879072746 

-0.8360956240 

0. 3879072746 

** 

17 

0.3879072746 

-0.8360956240 

-0.3879072746 

II 

18 

0.8360956240 

0.3879072746 

0.3879072746 

M 

19 

0.8360956240 

0.3879072746 

-0.3879072746 

" 

20 

0.8360956240 

-0.3879072746 

0.  3879072746 

" 

21 

0.8360956240 

-0. 3879072746 

-0.3879072746 

II 

B - 33.269905“ 


Table  2 - Direction  Cosines  and  Weights  for  2 x 25  Points. 


a 

•4 

*3 

ft 

w 

i 

1 

0 

0 

0.01269841058 

2 

0 

1 

0 

" 

3 

0 

0 

1 

II 

4 

0.7071067812 

0.7071067812 

0 

0.02257495612 

5 

0. 7071067812 

-0. 7071067812 

0 

II 

6 

0.7071067812 

0 

0. 7071067812 

(1 

7 

0.7071067812 

0 

-0.7071067812 

" 

8 

0 

0.7071067812 

0.7071067812 

»» 

9 

0 

0.7071067812 

-0.7071067812 

t» 

10 

0. 3015113354 

0. 3015113354 

0.9045340398 

0.02017333557 

11 

0.3015113354 

0.3015113354 

-0.9045340398 

• 1 

12 

0.3015113353 

-0.3015113354 

0.9045340398 

• 1 

13 

0.3015113354 

-0.3015113354 

-0.9045340398 

II 

14 

0.3015113354 

0.9045340398 

0.3015113354 

II 

15 

0.3015113354 

0.9045340398 

-0.3015113354 

II 

16 

0.3015113354 

-0.9045340398 

0.3015113354 

II 

17 

0.3015113354 

-0.9045340398 

-0.3015113354 

II 

18 

0.9045340398 

0.3015113354 

0.3015113354 

II 

19 

0.9045340398 

0.3015113354 

-0.3015113354 

II 

20 

0.9045340398 

-0.3015113354 

0.3015113354 

■ 1 

21 

0.9045340398 

-0,3015113354 

-0.3015113354 

ft 

22 

0.5773502692 

0.5773502692 

0.5773502692 

0.02109375117 

23 

0.5773502592 

0.5773502692 

-0.5773502692 

II 

24 

0.5773502692 

-0.5773502692 

0.5773502692 

II 

25 

0.5773502692 

-0.5773502692 

-0.5773502692 

II 

Pig.  3 

Comparison  with  some 
of  the  Test  Data  of 
Evans  and  Mara the 
(1968) 
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spacing  of  o-  directions  is  impossible. 

Hal-Ill l ami  Oil  ll?l  developed  numerical  Into- 
(ration  fenimlas  with  more  than  10  po into,  which 
give  consistent  results  even  in  the  Strain-soft- 
ening range.  The  most  efficient  formulas,  with  a 
nearly  uniform  spacing  of  a-d  fleet  ions,  are  ob- 
tained by  certain  subdivisions  of  the  sides  of  an 
Icosahedron  and/or  a dodecahedron  [17].  Much  for- 
mulas do  not  exhibit  orthoi;onul  symmetries. 

Other  formulas  wlhch  do  wore  also  developed  [17]. 
Taylor  series  expansions  on  a sphere  were  applied 
and  the  weights  associated  with  the  integration 

points  were  solved  from  the  condition  that  the 
greatest  possible  number  of  terms  of  the  Taylor 
series  expansion  would  cancel  out.  The  angular 
directions  of  certain  integration  points  were  fur- 
ther determined  so  as  to  minimize  the  error  term 
of  the  expansion,  formulas  involving  16,  21,  25. 

33,  37  and  61  points  were  derived,  with  errors  of 
8th,  10th  and  12th  order  [17).  Table  1 defines 
two  of  these  numerical  integration  formulas,  with 
21  and  25  points,  one  without,  and  one  with  ortho- 
gonal symmetry.  These  formulas  give  accuracy  that 
suffices  for  most  pratical  purposes,  for  crude 
calculations,  a formula  with  16  points  [17[  may 
sometimes  also  suffice.  The  directions  of  integra- 
tion points  are  illustrated  in  fig.  .’.  Also  shown 
one  stress-strain  diagrams  calculated  with  the 
formula  for  uniaxial  tension  in  various  directions 
with  regard  to  the  iutetgrat ion  points  (directions 
a,  b,  e,  d,,..)j  the  spread  of  the  curves  charac- 
terizes the  range  of  error. 


load  Increments,  and  subsequently  calculate  the  In- 
crements of  Cjj  and  o , . for  all  elements  and  all 

integration  points  in  each  element.  then  advance 
to  the  next  iteration  of  the  same  loading  stop, 
or  advance  to  the  next  loading  step- 


in  simulating  uniaxial  tensile  loading  of 
t ixed  direction,  the  unloading  cirterion  is  not 
important  since  the  only  unloading  occurs  at  mod- 
erate compressive  stresses,  for  which  a perfectly 
clastic  unloading  may  be  assumed. 

The  microplane  model  can  be  calibrated  by 
comparison  with  direct  tensile  tests  which  cover 
the.  strain-softening  response.  Such  tests,  which 
can  be  carried  out  in  a very  stiff  testing  machine 
and  on  sufficiently  small  test  specimens,  have 
been  performed  by  Evans  and  Marathe  [20]  ns  well 
as  others  [21-23].  Optimal  values  of  the 

three  parameters  of  the  model,  K , k,  and  n,  have 

been  found  [19]  so  as  to  achieve  the  best  fits  of 
the  data  of  Evans  and  Marathe.  Rome  of  these  fits 
are  shown  as  the  solid  lines  in  Eig.  3,  and  the 
data  are  shown  as  the  dashed  linos.  A better  test 
of  the  model  would,  of  course,  be  a tensile  test 
under  rotating  principal  stress  directions,  but 
such  tests  have  not  yet  been  performed. 

Note  that  with  this  theory,  one  has  only  two 
material  parameters,  E and  k,  to  determine  by 

fitting  test  data.  Trial  and  error  approach  is 
sufficient  for  that. 


Numerical  Algorithm 

The  following  numerical  algorithm  may  he  used 
for  the  mic.roplane  model  in  each  loading  step. 

1.  Determine  c^from  Kqs.  1 and  2 for  all 

n 

directions  a - l,...N.  In  the  first  iteration  of 
the  loading  step,  use  t ^ . for  the  end  of  the  pre- 
vious step,  and  in  subsequent  integrations  use  the 
value  of  tj.  determined  for  the  ir.id-stcp  in  the 

previous  iteration.  In  structural  analysis,  re- 
peat this  tor  all  finite  elements  and  for  ill  in- 
tegration points  within  each  finite  element. 

2.  For  ail  directions  evaluate  F1 2 3  (e  ) 

for  use  In  Kq.  7.  Also  check  for  each  direction 
whether  unloading  occurs,  as  indicated  by  viola- 
tion of  the  condition  s Ac  20.  If  violated, 
n n 

replace  F*  (e  ) with  the  unloading  stiffness  (which 
may  be  approximately  taken  as  K ^ ; however,  a 
better  expression  exists.) 

3.  Evaluate  D‘,  from  Eq.  7 and  D,  ,,  from 

ljkw  1 1 km 

Eq.  9.  In  structural  analysis,  repent  this  for 
all  elements  and  all  integration  points  In  each 
element.  When  solving  stress-strain  curves,  cal- 
culate then  the  increments  of  unknown  stresses 
and  unknown  strains  from  Eq.  9.  In  structural 
analysis,  solve  (by  the  finite  element  method)  the 
Increments  of  nodal  displacements  from  the  given 


Conclusion 

The  microplane  model  is  capable  of  simula- 
ting realistic  tensile  stress-strain  curve  with 
strain-softening  and  reduction  of  stresses  all  tlie 
way  to  zero.  Combined  with  the  blunt  crack  band 
concept,  In  which  the  strain-softening  is  re- 
stricted to  a region  of  a certain  characteristic 
width  that  is  a material  property  [l,2[,  this  model 
should  give  a realistic  representation  of  fracture. 
The  model  is  general  and  does  not  preclude  appli- 
cation to  stress  histories  in  which  the  principal 
stress  directions  rotate.  These  features  are 
particularly  attractive  for  the  analysis  of  tlm 
response  of  concrete  structures  subjected  to 
dynamic  loads. 
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